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Abstract 

We develop efficient group-theoretical approach to the problem of 
classification of evolution equations that admit non-local transforma- 
tion groups (quasi- local symmetries), i.e., groups involving integrals 
of the dependent variable. We classify realizations of two- and three- 
dimensional Lie algebras leading to equations admitting quasi-local 
symmetries. Finally, we generalize the approach in question for the 
case of an arbitrary system of evolution equations with two indepen- 
dent variables. 

1 Introduction. 

Consider the general evolution equation in one spatial dimension 

u t = F(t,x,u,u 1 ,u 2 , ■ ■ ■ ,u n ), n>2, (1) 

where u = u(t,x) is a real-valued function of two real variables t,x, Ui = 
d l u/dx l , i = 1, 2, ... ,n, and F is an arbitrary smooth real- valued function. 

The most general Lie transformation group leaving differential equation 
(CQ) invariant has the form (see, e.g., [HE]) 

t' = T(t,6), x' = X(t,x,u,6), u = U(t,x,u,6). (2) 
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Here T, X, U are smooth real- valued functions satisfying the non- singularity 
condition ^rff^P ^ in some open domain of M 3 and 9 = (9\, 62, ■ ■ ■ , 9 r ) £ 
W is the vector of group parameters. 

If a transformation of the space of variables t, x, u changes the specific 
form of ([1]) leaving invariant its differential structure, then we arrive at the 
concept of equivalence group. More precisely, if the (locally) invertible change 
of variables, 

t — > t' = T(t, x,u), x — > x' — X(t, x,u), u —> v! = U(t, x, u) 

maps Eq.([T]) into a possibly different n-th order evolution equation 

u' t , = G(t', x', u', u[, u' 2 , ...,u' n ), 

then this change of variables is called equivalence transformation. The set 
of all possible equivalence transformations forms a diffeomorphism group, S, 
and is called the equivalence group of Eq.([T]). Clearly, if we require that 
G = F then the equivalence group reduces to the symmetry group of Eq. (TTJ) . 
Consequently, Lie symmetry group of a given equation is a subgroup of its 
equivalence group. 

Let partial differential equation (JTJ) be invariant under Lie transforma- 
tion group ([2]). What would happen to this Lie symmetry if we perform a 
transformation from the equivalence group of the equation under study? Ev- 
idently, transformation group ([2]) after being rewritten in the 'new' variables 
t', x', vl becomes Lie symmetry of the transformed equation. 

Now suppose that we allow for the more general equivalence transforma- 
tion group 

t — > t' = T(t, x, u, v) , x — > x' = X(t, x, u, v), u — > u = U(t,x,u,v), 

where v = (u\, U2, ■ ■ ■ , u p , d~ l u, d~ 2 u, . . . d~ s u) with d^u = f u(t, x)dx and 
d^ k ^ 1 = d^ 1 d^ k . Saying it another way, we allow for an equivalence trans- 
formation to include derivatives and integrals of the dependent variable u. 
If such a transformation still preserves the differential structure of evolution 
equation (pQ), what would happen to Lie symmetries of the latter? The an- 
swer is, 'it depends'. In some cases, Lie symmetry transforms into another 
Lie symmetry. However, it could happen that some Lie symmetries would 
'disappear' after performing equivalence transformation, meaning that they 
cannot be found within the framework of the infinitesimal Lie method. The 
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reason is that the transformation rule for the variables t, x, u might con- 
tain derivatives and integrals of u, which are beyond reach of the standard 
Lie method. One needs to apply the generalized Lie [3]- [5] or non-Lie j6] 
approaches to be able to handle those symmetries. 

To single out Lie symmetries, which after non-local equivalence transfor- 
mation of an equation under study turn into non-Lie symmetries Ibragimov 
et al [7] introduced the term 'quasi-local symmetry', which we use through- 
out the paper. Independently, the concept of quasi-local symmetry has been 
suggested in [8]. 



In the present paper we suggest a simple regular method for deriving 
quasi-local symmetries (QLS) of evolution equations. Note that the basic 
idea of the method has been suggested in our paper [H] and some non-trivial 
examples of second-order evolution equations with QLS are given in [2] . 

Next, we demonstrate how to apply the method developed to arbitrary 
systems of evolution equations with two independent variables. 

2 Method description. 

The most general Lie transformation group admitted by evolution equation 
(CE]) is of the form (j2J). The infinitesimal operator, Q, of this group reads as 



that reduces Q to the canonical form d u (we drop the bars). Evolution 
equation (Tj[|) now becomes 



Note that the right-hand side of Eq.© does not depend explicitly on u. 
Differentiating (jSj) with respect to x yields 




(4) 



(3) 



u t = f(t,x,u 1 ,u 2 , • • .,«„)■ 



(5) 




n 



Making the change of variables 



(6) 
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and dropping the bars we finally get 



u -^l + ^l + Y^ df (7) 
dx du ^ dui-i 

where / = f(t,x,u,Ui, . . . ,u n _i). 

Thus non-local transformation ([6]) preserves the differential structure of 
the class of evolution equations (jSj). 

Let differential equation ((Sj) admit r-parameter Lie transformation group 
(121) with 9 = (6*1, . . . , 9 r ) and r > 2. To obtain the symmetry group of Eq.Q 
we need to transform (j2J) according to (jEJ). To this end we compute the first 
prolongation of formulas t^j and derive the transformation rule for the first 
derivative of u 

du' _ U u u x + U x 
dx' X u u x + X x ' 
So symmetry group (ED now reads as 

t' = T(t,6), x' = X(t,x,v,9), u' = " vU + U * (8) 

A. V U + Aj 

with v = d~ l u and U = (t,x,v,9). Consequently, if the right-hand sides of 
(JSJ) depend explicitly on the non-local variable v , then transformation group 
(JSJ) is a quasi- local symmetry of Eq.([7j). 

Evidently, transformations (jSj) include variable t> if and only if 

X v ^ or — — — — ^ 



<9v + 

or, equivalently (since all the functions involved are real- valued) , 



{X v ) 2 + (U VV X V — U V X VV ) 2 + {U XV X X — U X X X 



\2 



X t) J 



+ (U VV X X + U XV X V — U X X VV — U V X XV ) 7^ 0. (9) 

If X v ^ 0, then the above inequality holds true. If X v does vanish identi- 
cally, then (jHJ) reduces to U 2 V + U% v ^ 0. It is straightforward to express 
the above constraints in terms of the coefficients of the corresponding in- 
finitesimal operator of group ([2]). As a result, we get the following assertion. 
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Theorem 1 Equation (TJJ) can be reduced to evolution equation |?]) having 
QLS if it admits Lie symmetry, whose infinitesimal generator satisfies one 
of the inequalities 




du 

^ = 0, 
du 



(10) 




2 



) 
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^0. 



(11) 



Now we can formulate an algorithm for constructing evolution equations 
of the form ([1]) admitting QLS. 

1. We compute the maximal Lie symmetry group S of differential equation 



2. We classify inequivalent one-parameter subgroups of S and select sub- 
groups Si,...S p whose infinitesimal operators are of the form Q = 

£(t,x,u)d x + r](t,x,u)d u . 

3. For each subgroup, <Sj, we construct a change of variables (@J reducing 
the corresponding infinitesimal operator Q to the canonical operator 
du, which leads to evolution equations of the form (jSJ). 

4. Since the invariance group, S, admitted by (j3J) is isomorphic to S, 
we can utilize the results of subgroup classification of S. For each of 
the one-parameter subgroups of S we check whether its infinitesimal 
generators satisfies one of conditions ffTUl) . ffTTj) of Theorem 1. This 
yields the list of evolution equations that can be reduced to those having 



5. Performing non-local change of variables (El) yields evolution equations 
(J7j) admitting quasi-local symmetries (JSJ). 

Full implementation of the above approach will be the topic of our future 
publication. Here we restrict our considerations to classifying realizations 
of two- and three-parameter Lie transformation groups leading to evolution 
equations p]) that admit QLS. 

Hereafter we suggest that evolution equation ([1]) admits a Lie symmetry 
Q = £(t,x,u)d x + r](t,x,u)d u and therefore can be reduced to the form (JSJ). 
Differential equation (j5J) is guaranteed to admit at least a one-parameter 



QLS. 
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symmetry group, which is generated by the operator d u . What we are going 
to do is to describe all realizations of two- and three-dimensional Lie algebras, 
which 

• are Lie symmetry algebras of equations of the form (151), and, 

• have coefficients satisfying one of the inequalities ({TO]) . (flTI) from The- 
orem 1. 

With these realizations in hand, the problem of describing equations hav- 
ing QLS reduces to a straightforward application of the infinitesimal Lie 
method jU [TUJ [TT] , which boils down to integrating Euler-Lagrange system 
for calculating differential invariants of the corresponding Lie algebras of 
first-order differential operators. 

Let us remind that the most general symmetry generator admitted by 
([5]) is of the form ([3]) , while the most general equivalence group admitted by 
Eq.([T]) reads as (see, e.g., [16J) 

t = T(t), x = X(t,x,u), u = U(t,x,u), (12) 

where T, X, U are arbitrary smooth real- valued functions. 

We can always choose basis operators of Lie symmetry algebra of Eq. §5§ 
so that 

e = d u , ti = Ti(t)d t + £i(t, x, u)d x + r)i(t, x, u)d u , 

where i = 1, . . . ,r. Note that by the Magadeev theorem (12] the maximal 
possible value for r is n + 4, n being the order of evolution equation 
provided (JSJ) is not locally equivalent to a linear equation. In particular, for 
the second-order evolution equation we have r <= 6. By the definition of Lie 
algebra there are constant r x r matrix C and constant r-component vector 
c such that 

r 

[e , e»] = ^2 C iJ e 3 + c * e o> i = l,...,r. (13) 

3=1 

Here [Qi, Qz] = Q1Q2 - QiQ\- 

System of equations ffTB"]) is the starting point of our classification algo- 
rithm. First of all, let us note that by re-arranging the basis of the Lie 
algebra, e M — > X^=o c V e ^ Z 1 = 0, 1, . . . , r , we can always reduce the con- 
stant matrix C to the canonical form. Consequently, without any loss of 
generality we may assume that the matrix C is in the canonical form. 
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Computing commutators in the left-hand sides of (TT3|) and equating the 
coefficients of linearly-independent operators d t ,d x ,d u we get the following 
system of partial differential equations: 

«jL=Ofi l = Cn + c, Cf=0, (14) 

where £ = . . . , f r ), 17 = (in, . . . ,r) r ), f = (n, . . . , r r ). After integrating 
differential equations (|T4j) we need to ensure that the operators ei, . . . , e r do 
form a basis of Lie algebra and satisfy the additional set of commutation 
relations, 

r 

e j] = "^2 c ij e k, k = l,...,r. 

k=l 

Next, we simplify the form of operators e±, . . . , e r using the suitable equiv- 
alence transformations from the group S. As a final step, we verify that at 
least one of the coefficients of one the operators ei, . . . , e r satisfy either (flUi) 
or (TUD. 



3 Realizations of two-dimensional QLS alge- 
bras 

For the case when r = 1, system (j!4p reduces to a pair of non-coupled 
differential equations 

£u = Af , r] u = Xr] + c, At = 0, (15) 

where A, c are constants. 

While integrating ffl5l) we need to differentiate between the two cases 
A ^ and A = 0. 

Case 1. A 7^ 0. The general solution of (|15p has the form 

r(t) = 0, £(t,x, u) = Wi(t,x) exp(Aw), 
r)(t, x, u) = W 2 (t, x) exp(Aw) — cA _1 . 

where Wi,W2 are arbitrary smooth functions. So that the two-dimensional 
Lie algebra (e , ex) read as 

Co = d u , ei = W±(t, x) exp(Xu)d x + [W2(t, x) exp(Aw) — cA^ 1 ) d u . 
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As A 7^ we can always re-scale u, i.e., make a transformation u — *> ku, k = 
const, in order to get A = 1. Next taking as new e\ the linear combination 
cA -1 eo + e\ we get rid of the term in e± which is proportional to c, namely, 

e = d u , ei = Wi(t, x) exp(u)d x + W 2 {t, x) exp(u)d u . 

It is not difficult to verify that the most general subgroup of the equivalence 
group S not altering the form of equation (jSJ) and operator d u is given by the 
formulas 

t = T(t), x = X(t,x), u = u + U(t,x), (16) 

where T, X, U are arbitrary smooth functions. 

Since the functions W± and W 2 do not vanish simultaneously, we have 
three possible subcases, (1) W x ^ 0, W 2 ^ 0, (2) W x ^ 0, W 2 = 0, and 
Wi = 0,W 2 ^ 0. 

Case 1.1. W\ 7^ 0, W 2 ^ 0. Applying (Tl6l) with T = t we reduce the operator 
e\ to the form 

e x = ei expfw)^ + e 2 exp(w)<9 u , 

where e\ = ±1,62 = ±1- Combining equivalence transformation £ — > t, a; — > 
— x, w — > w and re-scaling ei — > — ei we get the final form of the basis 
elements eo, ei 

e = <9 U , ei = exp(u){d x + d u ). 

Since £ u = exp(u) ^ 0, the condition (jTDJ) of Theorem 1 holds true and 
the evolution equation invariant under the above algebra is equivalent to a 
quasi-linear evolution equation that admits QLS. 

Case 1.2. Wi ^ 0,W 2 ^ 0. Applying transformation flT3D with T(£) = 
t,U(t,x) = 0, = 1 we reduce the operator ei to the form ei = exp(u)d x . 
Again, the coefficient £ = exp(w) obeys condition (TTOT) of Theorem 1 and, 
therefore, it gives rise to the two-dimensional Lie algebra 

e = d u , ei = exp(u)d x 

that leads to an evolution equation admitting QLS. 

Case 1.3. W\ — 0, W 2 ^ 0. Applying transformation fTToD with T(t) = 
t,X(t,x) = x we reduce the operator e± to the form ei = ± exp(w)<9 u . Re- 
scaling, if necessary, the operator e± to —e\ we can make sure that e± reads 
as exp(w)<9 u and finally get 

e = d u , e 2 = exp(u)d u . 
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Since the coefficient r\ = exp(w) obeys condition fill I) of Theorem 1, an evo- 
lution equation invariant under the above algebra is equivalent to a partial 
differential equation of the form (|7]) which has QLS. 
Case 2. A = 0. System (fT5l) is readily integrated to yield 

T(t) = W Q (t), £ = W 1 (t,x), r) = W 2 (t,x) + cu. 

Here Wq, Wi, W 2 are arbitrary smooth functions. Checking the conditions of 
Theorem 1 we see that neither of them can be satisfied by the coefficients of 
the operator e x . Consequently, this case yields no equations admitting QLS. 

Below we give the full list of £-inequivalent realizations of two-dimensional 
Lie algebras spanned by the operators eo = d u , e± = T{t)d t +£(t,x,u)d x 
+rj(t,x,u)d u satisfying the conditions of Theorem 1 

A\ : (d u , exp(u)(d x + d u )), 
A\ : (d u , exp(u)d x ), 
A\ : (d u , exp(u)d u ). 

Evolution equations ([5]) invariant under the above algebras are reduced to 
differential equations that admit QES. The corresponding QLS are obtained 
by re-writing the transformation groups generated by the operators e\ in 
terms of new (non-local) variables t, x, u and d^u. 

Consider, for example, the algebra A\ = (d u , exp(u)d x ). Applying the 
standard infinitesimal Lie algorithm [10] we obtain the determining equations 
for the function / 

-«J - fx + u 2 J Ux + {u 3 x + 3u x u xx )f Uxx = 0. 

The general solution of the above equation reads as 

fit, x, u x , u xx ) = U x f((J , UJi, u 2 ), 

where / is an arbitrary smooth function and u> = t, uj\ = (xu x — l)^" 1 , u> 2 = 
[u xx + u^)u~ 3 are absolute invariants of the transformation group generated 
by the operators d u and exp(u)d x . Consequently, the evolution equation 
invariant under the algebra A\ is of the form 

Ut = Uxf (^0,^1,^2) ■ 
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Differentiating the above equation with respect to x and replacing u x with u 
according to we arrive at the evolution equation 

~ u x + u 2 ~ uu xx - 3(w 2 + l)u x ~ 

u t — u xJ H 5 — Jw! H j Jw 2 

u z u 

with ujq = t, lo\ = x — tt _1 , u 2 = {u x + u 2 )u~ 3 . This differential equation 
admits the following quasi-local symmetry group 

u 

t' = t, x' = x + 8 exp(v) , u 



1 + 6u exp(w) 
where 6 is a group parameter and v = d^u. 



4 Realizations of three-dimensional QLS al- 
gebras 

Consider system of partial differential equations (1141) with r = 2. The con- 
stant 2x2 matrix C has been reduced to the canonical real Jordan form. 
There are three inequivalent cases that need to be considered separately, 
namely, when eigenvalues 

1. are complex conjugate, 

2. are real, and the matrix C is diagonal matrix, 

3. are real and the matrix C is the 2x2 canonical Jordan box 




We consider in detail the class of realizations of three-dimensional Lie alge- 
bras obtained for the case when C has two complex eigenvalues Ai, A2. For 
the remaining two classes we present the final results only. 



4.1 Case of diagonal canonical form with complex ei- 
genvalues 

As the characteristic equation of the real matrix C is real, the eigenvalues 
have to satisfy the additional constraint A* = A2- Consequently, if we define 
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A = (Ai + A 2 )/2 and a = (Ai — A 2 )/(2i), then the general solution of <HM can 
be represented in the form 

n = 0, t 2 = o, 

£1 = (Wi(t, x) cos(cra) + W 2 (t, x) sin(cra)) exp(Aw), 

£2 = (W 2 (t, x) cos(a«) — x) sin(ora)) exp(Au), (17) 

771 = (W 3 (t, x) cos(au) + W^t, x) sin(au)) exp(Aw) + 61, 

rj2 = (W^(t, x) cos(aw) — Wz(t, x) sin(att)) exp(Aw) + 6 2 - 

Here W\, W 2 , W3, W4 are arbitrary smooth real- valued functions, 61,62 are 
real constants. Hence, the most general form of the basis operators e±, e 2 is 

t\ = (Wi(t, x) cos(au) + W2(t, x) sin(em)) exp(Xu)d x 

+ (^{W 3 (t,x) cos(aw) + W±(t,x) sin(ait)) exp(Aw) + bAd u , 

e 2 = (W2(t, x) cos(aw) — Wi(t, x) sin(aw)) exj>(\u)d x 

+ ((W^t, x) cos(au) — Ws(tj x) sin(em)) exp(Aw) + 62 \d u . 

Note that a 7^ 0, otherwise, Ai, A2 are not complex. By re-scaling the variable 
u we can make a equal to 1. Next, applying to the operators ei,e 2 an 
equivalence transformation (|15p with T(t) = t,X(t,x) = x we can get rid of 
the function W 2 . With these remarks the operators e±, e 2 take the form 

e\ = ((W3(t, x) cosu + Wi(tj x) sinw) exp(Aw) + b\\d u 

+Wi(t, x) cos uexp(Xu)d x , (18) 

e 2 = (^(Wi(t,x) cos u — Ws(t,x) sin u) exp(Xu) +b2jd u 
—Wi(t, x) sin-u exp(Xu)d x . 

In what follows we need to distinguish between the cases, A 7^ and A = 0. 
Case 1. A / 0. Performing, if necessary, transformation (|15p with T(t) = 
t, U(t, x) = u we can always make non-vanishing identically function W equal 
to 1. Next, taking as e% and e 2 the linear combinations e% — b\eo and e2 — 62^0 
we can get rid of parameters 61, 62- 

Now we need to ensure that the operators e , e 1; e 2 do form a realization 
of a Lie algebra. To this end we have to verify that the relation 

[ei, e 2 ] = ae 1 + (3e 2 + 7e (19) 
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with some real a, /3, 7 holds true. Calculating the commutators and equating 
the coefficients of linearly-independent operators d t ,d x , d u we get the system 
of differential equations for W3, W4. Its general solution is given by the 
formulas 

W 3 = A(A 2 + l)" 1 ^ + pit))' 1 , W 4 = -(A 2 + l)" 1 ^ + p(t))-\ 

Here p(t) is an arbitrary smooth function. Making the equivalence trans- 
formation (|15p with T(t) = t,X(t,x) = p(t),U(t,x) = we eliminate the 
function p(t) and arrive at the following realization of a three-dimensional 
Lie algebra 

e = d u , 

e% = exp(Aw) cosud x + (A 2 + 1) _1 (A cosw — sin-u)x _1 exp(Aw)<9 u , 
e 2 = — exp(Aw) siaud x — (A 2 + 1) _1 (cosm — A sin^x^ 1 exp(Aii)<9 n , 

Evidently, the coefficients of ei,e 2 satisfy condition fflOl) of Theorem 1 and, 
consequently, evolution equation invariant under the symmetry algebra eo, 
ei, e 2 can be reduced to the one having QLS. 
Case 2. A = 0. Operators ffTS]) take the form 

e\ = ({W3 {t, x) cos u + W^i, x) sin u) + 61^ d u + W\ (t, x) cos ud x , 
e 2 = M Wi(t, x) cosu — W%(t, x) sinw) + b 2 jd u — Wi(t, x) sinud x . 

Replacing e\ and e 2 with the linear combinations e% — bieo and e 2 — 6 2 eo 
eliminates the parameters 61,62- So that we can assume that hi = 0, 6 2 = 
without any loss of generality. 

Case 2.1. W\ 7^ 0. Utilizing equivalence transformation ffl5l) with T = 
t, U = we can make W\ equal to 1. After simple algebra we prove that for 
the operators ei,e 2 to satisfy the remaining commutation relation flTTTT) the 
functions W3, W4 have to take one of the following forms: 

W 3 = 0, W 4 = /itan(/ix), 
Wz = 0, W4 = — /itanh(/ix), 
W 3 = 0, W 4 = x~\ 

where fi is an arbitrary real parameter. Inserting the above expressions into 
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the corresponding formulas for ei, e 2 we finally get 



e = d u , 

e\ = cos ud x + fi ta>n(iAx) sin ud u , 

e 2 = — sin ud x + /x tan(/ix) cos ud u ; 

e = 9«, 

ei = cos ud x — /itanh(/ix) sin-u<9 n , 

e 2 = — sin wc^. — /i tanh(/ix) cos w<9 n ; 

e = d u , 

e\ = cosud x — x~ l sinud u , 

e 2 = — sinud x — x~ x cosud u . 

Case 2.2. W\ = 0. In this case using transformation (|15|) with T = t, X = x 
we can eliminate W4. Inserting the corresponding expressions for ei,e2 into 
(|T9|) and solving the obtained equations within the equivalence relation S 
yields 

Co = d u , t\ = cos ud u , e 2 = smud u . 

Note that all realizations of three-dimensional Lie algebras obtained under 
Cases 2.1, 2.2 satisfy the conditions of Theorem 1. Consequently, evolution 
equations invariant with respect to the above algebras can be transformed 
into equations admitting QLS. 

Summing up we present the full list of realizations of three-dimensional 
Lie algebras, obtained for the case when 2x2 matrix C in (|14p has two 
complex eigenvalues. 

A\ : (d u , exp(/iw) cosud x + (/i 2 + 1)~ (/icosu — sinu)x^ exp(fiu)d u , 
— exp(/iw) smud x — (/i 2 + 1)~ (cosu — /isin^x -1 exp(/xtt)<9 u ), 
Al : (d u , cosud x + /xtan(/xx) smud u , — smud x + /xtan(/xx) cosud u ), 
A\ : (d u , cosudx — /itanh(zxr) smud u , — sinud x — /itanh(/ix) cosud u ), 
A\ : (d u , cosudx — x" 1 smud u , — sin ud x — x~ x cos ud u ) , 
Al : (d u , cosud u , sin ud u ). 

Here \x is an arbitrary real constant. 

Evolution equations jSJ) invariant under the algebras A3, . . . , A\ can be 
reduced to equations admitting QLS. 
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4.2 Case of diagonal canonical form with real eigenval- 
ues 

Without loss of generality we may assume that the matrix C from (|14[) has 

/ Ai \ 

been reduced to the diagonal matrix I . Since Ai, A 2 do not vanish 

\ A 2 J 

simultaneously we may assume that Ai 7^ 0. Re- scaling u we make Ai equal 
to 1. With this choice of C, system ( 114)) takes the form 

7i = 0, A 2 T 2 = 0, 

£iu = £,1, £,2u = A 2 £ 2 , 

Viv = Vi + ci, V2u = A 2 ?7 2 + c 2 . 

Integrating the above system within the equivalence relation £ , inserting the 
result into the remaining commutation relation (Tl9l) and solving the equa- 
tions obtained we arrive at the following realizations of three-dimensional Lie 
algebras: 





(d u , 


exp(u)d x , (x 2 + cr(t)) exp(— u)d x + 2xexp(—u)d u ), 


Al 


(d u , 


exp(u)d x + eexp(u)d u , (cr(t) exp(x) ± exp(2x) + /1) exp(— u)d x 




+ (±exp(2x) - /i)d u ), 


AI 


(d u , 


exp(u)d x + exp(u)d u , (cr(£) exp(— fix) ± exp((l — [x)x)) 




X 


exp(fiu)d x ± exp((l - /i)x)) exp(/iu)d u ), 


Al 


(d u , 


exp(u)d x , a(t)x exp([m)d x + a{t) exp(fm)d u ) , 


Al° 


(d u , 


exp(u)d x , a(t)exp([m)d u ), 


Al 1 


(d u , 


exp(u)(d x + d u ), {cr{t) + exp(x))d x + exp(x)d u + ed t ), 


A? 


(d u , 


exp(u)(d x + d u ), cr(t)d x + ed t ), 


Af 


(d u , 


exp(u)d x , a(t)(xd x + d u ) + ed t ), 


Af 


(d u , 


exp(u)d x , a(t)(d x + d u ) + ed t ), 


Af 


(d u , 


exp(u)d u , cr(t)d u + ed t ), 


Af 


(d u , 


exp(u)d u , d x + cr(t)d u ), 


Al 7 


(d u , 


exp(u)(d x + d u ), exp(u)(a 1 (t) exp(-x) + a 2 {t))d x 




+ 


01 (t) exp(u)<9„) , 


Al 8 




exp(u)d x , exp(M)(cri(t)x + a 2 {t))d x + ai(t) exp(u)d u ), 
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Af : (d u , exp(u)d u , xexp(u)d u ), 
Af : (d u , exp(u)d u , texp(u)d u ). 

Here u(t), cxi(t), 02 (i) are arbitrary real- valued smooth functions, // is an 
arbitrary real parameter, and e = 0, 1. 

By construction the coefficients of algebras A3 — A3 satisfy the conditions 
of Theorem 1. Consequently, evolution equations (jSJ) invariant under these 
algebras can be reduced to ones having QLS. 

4.3 Case of non-diagonal canonical form 

For the case under consideration, the matrix C from fTH|) is of the form 




. System of partial differential equations (jl"3l) now reads as 



ATi + T 2 = 0, XT 2 = 0, 

?7i« = Ar/i + 772 + ci, r]2u = Ar] 2 + c 2 . 



Integrating the above system within the equivalence relation £, inserting the 
result into (1191) and solving the relations obtained we obtain eleven realiza- 
tions of three-dimensional Lie algebras: 



Af 


(d u , 


uexp(u)d x + x _1 (-u + 1) exp(u)d u , exp(u)d x + x" 1 exp(w)<9 u ) 


Af 


(d u , 


ud x + (fiu 2 + a (t) exp(4fix))d u , d x + 2fid u ), 


Af 


(d u , 


ud x + (fix + vu + oi(t))d u + a 2 (t)d t , d x ), 


Af 


(du, 


d x + fiu tan(/ix)<9 u , tsn(fix)d u ) , 


Af 


(du, 


d x — fiuta.nh(fix)d u , tanh(/ix)<9 u ), 


Af 


(d u , 


d x x zid u , x d u ) , 


Af 


(du, 


(w 2 + x)d u , ud u ), 


Af 


(d u , 


(u 2 + t)d u , ud u ), 


Af 


(d u , 


u[y + 2/xtan(/xt + ai))9 n + 2d t , [y + 2/itan(/it + ax))<9 u ), 


Af 


(du, 


w(z/ + 2/itanh(— fit + ax))<9 u + 2<9 i; [y + 2/i 




x tanh(— fit + ax))d u ), 


Af 


(d u , 


(/iz/x — fit — 2){t — vx) 1 ud u + 2td t , (fivx — fit — 2) 




X 


(t — h>x)~ l ud u ). 
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Here cr(t), <Ji(t), <y 2 {t) are arbitrary real- valued smooth functions, a,fi,i> are 
arbitrary real parameters. 

Evolution equations (j3J) admitting symmetry algebras A% 1 , . . . , Af 1 can 
be reduced to equations having QLS. 

5 Some generalizations 

The technique developed in the previous sections naturally expands to cover 
general systems of evolution equations 

u t = F(t,x,u,ui, . . . ,u n ), n>2. (20) 

Here F = (F 1 , . . . , F m ) is an arbitrary m-component real-valued smooth 
function, u = (u 1 , . . . , u m ) G M m , and Ui = (d l u) / (dx l ), i — 1, . . . , n. 

Suppose that system of evolution equations (1201) admits m-parameter 
Abelian symmetry group which leaves the temporal variable, t, invariant. 
The infinitesimal generators of this group have to be of the form 

m 

Qi = x, u)d x +y] ijjj {t, x, •u)9 tt j , i = l,...,n (21) 

3=1 

and what is more, the rank of the matrix composed of the coefficients of 
differential operators d x , d u i, . . . , d u ™ equals to m. Given these conditions, 
there exists a change of variables 

i=t, x = X(t,x,u), u = U(t,x,u) (22) 

that reduces operators ff2T|) to canonical ones Qi = d U i, i = 1, . . . ,m (see, 
e.g., [ID])- Consequently system of evolution equations (12171) takes the form 

u t = f{t,x,u u . . . ,u n ), n>2. (23) 

Note that we dropped the bars. 

Now we can apply the same trick we utilized for the case of a single 
evolution equation. Namely, we differentiate (1231) with respect to x and map 
u — > u x thus getting 

*=s + EE^ (24) 

i=i j=i *— j- 
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with u 3 Q = uK 

Now if the system of evolution equations (1231) admits the Lie transforma- 
tion group 

t' = T(t,9), x' = X{t,x,u,6), u' = U(t,x,u,9), (25) 

where 9 G IR is a group parameter, then the transformed system of equations 
(I2"5"j) admits the group 

t> =T(t,9), 

x' = X(t,x,v,9), (26) 

with v l = d~ l u l = f u % dx and U = U(t,x,v,9). Consequently, provided 
either of relations 

holds for some j, 1 < j < m, system of evolution equations admits QLS 

(USD. 

Set of relations (|2"T|) is equivalent to the following system of inequalities 

m 

E x l + °> 

i=l 

or 

m mm 

E^ = ' £(^) 2 + E (4)Vo. 

i=l *J=1 i,j,k=l 

Note that all the functions involved are real- valued, so that vanishing of the 
sum of squares requires that every summand vanishes individually. Rewriting 
the obtained relations in terms of coefficients of the corresponding infinites- 
imal operators we arrive at the following assertion. 

Theorem 2 System of evolution equations (TJSp can be reduced to a system 
having QLS if it admits Lie symmetry whose infinitesimal operator Q = 
r(t)d t + £(£, x, u)d x + Y^iL\ Vii^, x, u)d u i satisfies one of the inequalities 

m 

E & / °> (28) 
i=i 
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7)1 



rn 



in 



(29) 



Summing up we formulate the algorithm for classifying systems of evolu- 
tion equations fl20|) admitting QLS. 

1. We compute the maximal Lie symmetry group S of system of partial 
differential equations (1201) . 

2. We classify inequivalent m-parameter Abelian subgroups Si,...S p of 
the group S and select subgroups whose infinitesimal operators are of 
the form fT2T]) . 

3. For each subgroup Si we construct change of variables (1221) reduc- 
ing commuting infinitesimal operators, Qi, to the canonical forms <9 S i, 
which leads to system of evolution equations (1231) . 

4. Since the invariance group, S, admitted by ( 1231) is isomorphic to S, we 
can utilize the results of subgroup classification of S. For each of the m- 
parameter Abelian subgroups of S we check whether their infinitesimal 
generators satisfy one of conditions (1281) . (129]) of Theorem 2. This yields 
the list of systems of evolution equations that can be reduced to those 
having QLS. 

5. Performing the non-local change of variables u l — > u l x , i = l,...,m 
we obtain systems of evolution equations ( 1241) admitting quasi-local 
symmetries ( |26l ). 

We intend to devote a special publication to application of this algorithm 
to Schrodinger-type systems of partial differential equations. Here we present 
an example of Galilei-invariant nonlinear Schrodinger equation, which leads 
to the equation possessing QLS. 

Consider nonlinear Schrodinger equation 



where V = i>RE(t,x) + iipi M {t,x), V = ^RE(t,x) - iip IM (t,x), a ^ is 
an arbitrary real constant and F is an arbitrary complex-valued function. 



#t = V'axc + + i<x) Vx - (V 2 )(a; + ia) 
+F(2ia(x + ia)ip x — (x — ia)(ip — ip*)), 



(30) 
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According to [13] , Eq. (130 j) admits the Lie algebra of the Galilei group having 
the following basis operators: 

ei = <9 t , 

e 2 = 9^, + <V , 

e 3 = (x + ia)~ l d^ + (x — ia) , 

e4 = — (t + (x + ia)~ l il))d^, — (t + (x — ict) ij)*)d^* . 

Operators e 2 , e% commute and the rank of the matrix of coefficients of 
operators d t ,d x ,d,p,d^* is equal to 2. Consequently, there is a change of 
variables that reduces e 2 , e% to canonical forms d u , d v . Indeed, making the 
change of variables 

u{t,x) = (l/2)(ip+ip*), v(t,x) = (2ia)- 1 (x 2 + a 2 )(ij-tp*) (31) 

transforms e 1; e 2 to become t\ = d u , e 2 = d v . So we can apply the above 
algorithm to Eq. (l3"U|) transformed according to fl3"T|) . As the coefficients of 
the transformed operator e 3 satisfy f[2"9~j) . it leads to QLS of the system of 
evolution equations of the form fl2M . 

6 Concluding remarks 

In the present paper we develop the efficient approach to constructing evolu- 
tion type partial differential equations which admit quasi-local symmetries. 
It is important to emphasize that the algorithm can be applied iteratively. 
Namely, if the transformed equation possesses Lie symmetries which satisfy 
conditions of Theorem 1, then it again can be transformed to a new evolution 
equation admitting QLS and so on. What is more, the equation obtained 
as the iVtli iteration of the algorithm admits QLS which involves non-local 
variables <9 _1 w, d~ 2 u, . . . , d~ N u. 

It is also of great interest to explore the case of multi-component evolution 
equations. The most natural objects are the non-linear Schrodinger-type 
equations or systems of nonlinear reaction-diffusion equations. 

There is a different approach to analyzing non-local symmetries for some 
special differential equations based on the notion of potential symmetries 
introduced by Bluman O [15]. It is of interest to explore the connection 
between two approaches, which would open a way to group classification of 
equations with potential symmetries. 
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Study of the above mentioned problems is in progress now and will be 
reported in our future publications. 
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